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ABSTRACT 

We acquired near-infrared spectra of the Oort cloud comet C/2012 F6 (Lemmon) at three different heliocentric 
distances ( R h) during the comet’s 2013 perihelion passage, providing a comprehensive measure of the outgassing 
behavior of parent volatiles and cosmogonic indicators. Our observations were performed pre-perihelion at 7? h = 

1.2 AU with CRIRES (on 2013 February 2 and 4), and post-perihelion at R h = 0.75 AU with CSHELL (on March 
31 and April 1) and R h = 1.74 AU with NIRSPEC (on June 20). We detected 10 volatile species (H 2 0, OH* prompt 
emission, C 2 H 6 , CH 3 OH, H 2 CO, HCN, CO, CH 4 , NH 3 , and NH 2 ), and obtained upper limits for two others (C 2 H 2 
and HDO). One-dimensional spatial profiles displayed different distributions for some volatiles, confirming either 
the existence of polar and apolar ices, or of chemically distinct active vents in the nucleus. The ortho-para ratio for 
water was 3.31 zb 0.33 (weighted mean of CRIRES and NIRSPEC results), implying a spin temperature >37 K at 
the 95% confidence limit. Our (3cr) upper limit for HDO corresponds to D/H < 2.45 x 10 3 (i.e., < 16 Vienna 
Standard Mean Ocean Water, VSMOW). At R^ = 1.2 AU (CRIRES), the production rate for water was 2(H 2 0) = 

1.9 =b 0.1 x 10 29 s _1 and its rotational temperature was T mi ~ 69 K. At R h = 0.75 AU (CSHELL), we measured 
2(H 2 0) = 4.6 zb 0.6 x 10 29 s- 1 and T rot = 80 K on March 31, and 6.6 zb 0.9 x 10 29 s- 1 and T Iot = 100 K on April 1. 

At R h = 1.74 AU (NIRSPEC), we obtained g(H 2 0) = 1.1 zb 0.1 x 10 29 s- 1 and T mi ~ 50 K. The measured volatile 
abundance ratios classify comet C/2012 F6 as rather depleted in C 2 H6 and CH 3 OH, while HCN, CH 4 , and CO 
displayed abundances close to their median values found among comets. H 2 CO was the only volatile showing a 
relative enhancement. The relative paucity of C 2 H6 and CH 3 OH (with respect to H 2 0) suggests formation within 
warm regions of the nebula. However, the normal abundance of HCN and hypervolatiles CH 4 and CO, and the 
enhancement of H 2 CO, may indicate a possible heterogeneous nucleus of comet C/2012 F6 (Lemmon), possibly 
as a result of radial mixing within the protoplanetary disk. 

Key words: astrochemistry - comets: general - comets: individual (C/2012 F6 (Lemmon)) - molecular 
processes - Oort Cloud - planets and satellites: formation 
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1. INTRODUCTION 

Comets visit the inner solar system from at least two long- 
term reservoirs — the scattered disk/Kuiper Belt (KB) and the 
Oort cloud (OC). But they actually agglomerated in other re- 
gions during the formative phase of the solar system. Indeed, 
icy bodies (such as comets and trans-Neptunian objects) are be- 
lieved to preserve material from that time. While some may have 
experienced subsequent processing owing to the effects of ther- 
mal warming, exposure to energetic radiation, and/or collisions 
since their formation, the relatively pristine composition of oth- 
ers may reveal the physicochemical conditions that prevailed 
4.6 billion years ago. 

Although an ongoing subject of discussion, current evidence 
suggests that most comets might have formed in an extensive 
region near or around the C0 2 and CO snow lines (over the range 
in heliocentric distance R h ~ 5-30 AU; Dodson-Robinson et al. 
2009; A’Hearn et al. 2012; Qi et al. 2013) before migrating to 


* Based on observations obtained at the European Southern Observatory at 
Cerro Paranal, Chile, under program 290.C-5016, at the NASA’s Infrared 
Telescope Facility, program 2013A_071, and Keck Observatory, program 
H233NS, at Mauna Kea, Hawaii. 


their current reservoirs (OC and KB; e.g., Gomes et al. 2005; 
Brasser & Morbidelli 2013). Further lines of evidence also 
suggest that a stationary formation (i.e., one not accounting 
for radial mixing) might not fully account for the diversity 
of volatile abundances and other cosmogonic parameters in 
icy bodies, such as spin temperatures and isotopic ratios (see 
Paganini et al. 2012, 2013a, and references therein). Indeed, 
prior to outward migration, small bodies could have experienced 
possible chemical processing as a result of dynamical translation 
within the disk driven by the inward and outward migration of 
Jupiter and Saturn in the first 100,000 yr after formation of the 
Sun (the so-called Grand Tack event; cf. Walsh et al. 2011 and 
references therein). 

Thus far, the analysis of primary volatiles in cometary comae, 
i.e., species directly released from the nucleus (traditionally 
termed parent species ), represents unique evidence for (rela- 
tively) unaltered material. The importance of such measure- 
ments is twofold: it might help us unveil the primordial condi- 
tions that occurred during the early solar system, and it might 
help us prove whether dynamical conditions played an important 
role during planetary formation. In addition to studies of product 
species at optical wavelengths (e.g., A’Heam et al. 1995; Fink 
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Table 1 

Log of the Observations 


Paganini et al. 


Instrument 

Date 

(UT 2013) 

Rh 

(AU) 

A 

(AU) 

A-dot 
(km s -1 ) 

P.A. 

(°) 

a 

(°) 

Settings 

(Time on Source [minutes]) 

CRIRES 

Feb. 2 

07:28-09:32 

1.20 

0.99 

-4 

262.4 

53.4 

HCN(32), CH 3 OH(12), CH(20) 

CRIRES 

Feb. 4 

05:25-06:40 

1.18 

0.99 

-2 

287.6 

52.3 

H 2 O(20), HDO(24) 

CSHELL 

Mar. 31 (Apr. 1) 
19:45-(00:28) 

0.75 

1.55 

+ 16 

199.3 

32.4 

CO(12), C 2 H 6 (12), H 2 0(32), CH 4 (40), H 2 CO(32) 

CSHELL 

Apr. l(Apr. 2) 
19:46-(01:08) 

0.75 

1.56 

+ 15 

201.8 

31.8 

C 2 H 6 (16), CH 3 OH(26), H 2 0(32), CO(12) 

NIRSPEC 

Jun. 20 
12:04-15:18 

1.74 

1.79 

+5 

256.3 

33.4 

KL1(40), KL2(76) a 


Notes. R h: heliocentric distance; A: geocentric distance; A-dot: geocentric velocity; P.A.: position angle of the extended Sun-comet vector; or. solar 
phase (Sun-comet-Earth) angle. (These values represent the midpoint of data acquisition.) 

a Both in KL1 and KL2 we targeted H 2 O, C 2 H 6 , and CH 3 OH, while KL2 also allowed estimates of H 2 CO, HCN, NH 3 , NH 2 , and C 2 H 2 . 


2009; Schleicher & Bair 2010; Langland-Shula & Smith 2011; 
Cochran et al. 2012), radio and infrared surveys of primary 
volatiles seek to build taxonomical databases for comets (for 
reviews, see Bockelee-Morvan et al. 2004; Mumma & Chamley 
2011). Most of these studies show compositional diversity 
amongst objects from the OC and KB, thus demonstrating the 
need for more comprehensive statistical evaluations. Aided by 
dynamical (theoretical) studies and sample analyses of solid ma- 
terial (e.g., of returned dust samples; Brownlee et al. 2006), our 
ongoing IR survey is providing clues into the possible origins 
of the early solar system by analyzing a growing database of 
comets from both the OC and KB. 

On 2012 March 23, Alex Gibbs of the Mount Lemmon Sur- 
vey discovered an inbound faint object with a visual magnitude 
of ~20 at 5 AU from the Sun, which subsequently was desig- 
nated C/2012 F6 (Lemmon; hereafter C/2012 F6). Its eccen- 
tricity ( e = 0.998), original semimajor axis 6 ( a = 458.08 AU), 
Tisserand parameter (7} = 0.147), and inclination to the eclip- 
tic (i = 82? 6) classify it as belonging to the nearly isotropic 
(long-period) comet group originating from the OC (Nakano 
2013). In January 2013, the comet underwent significant bright- 
ening. This prompted several optical observations, which re- 
vealed a greenish coma, indicative of a volatile-dominated 
spectrum (e.g., atomic oxygen and various carbon- and 
nitrogen-bearing species). We initiated a worldwide campaign 
using different facilities, including VLT/CRIRES, Keck-2/ 
NIRSPEC, and IRTF/CSHELL, in the infrared. Here, we 
present results obtained for several volatiles in C/2012 F6 dur- 
ing this campaign, including production rates, rotational temper- 
atures, spatial profiles, and cosmogonic parameters (ortho-para 
ratio and D/H in water). 

2. OBSERVATIONS AND DATA ANALYSIS 

We obtained high-resolution near-IR spectra of comet 
C/2012 F6 on five dates spanning February through June 2013 
(bracketing perihelion, q = 0.73 AU on 2013 March 24), using 
three separate spectrometer-telescope combinations. Table 1 
presents a log of our observations. 

With CRIRES we used the 0?4 slit, which delivers high spec- 
tral resolving power (X/AX ~ 50,000) and spatial coverage of 


6 S. Nakano Note, NK 2452, http://www.oaa.gr.jp/~oaacs/nk/nk2452.htm 


40". The adaptive optics module (Multi- Applications Curva- 
ture Adaptive Optics) provided seeing correction and image 
stabilization, thereby delivering excellent spatial registration. 
On both nights, we oriented the slit at position angle 290°, i.e., 
close to the Sun-comet direction (262° on February 2 and 288° 
on February 4). Atmospheric conditions were stable on both 
nights: relative humidity was ~20%, wind speeds were below 
7 ms -1 , and seeing was in the range of 1"0-1"5. Standard star 
HR 5336 (located near the comet’s trajectory) allowed flux cal- 
ibration and measures of column content of absorbing species 
in the terrestrial atmosphere. 

With CSHELL we used the 2 " wide slit, maximizing the 
flux collection while still delivering moderate resolving power 
(X/AX ~ 20,000), which was adequate for our molecular 
studies. This slit width was necessary because C/2012 F6 
was positioned 20°-25° south of the Sun and so was available 
only during daytime, which precluded use of the CCD guider 
in CSHELL. The 30" long slit was oriented north-south on 
March 31, close to the projected Sun-comet direction (position 
angle 199°; Table 1), and east-west on April 1. Flux levels the 
terrestrial contribution were established through observations 
of the standard star HR 8728. For our CSHELL observations, 
three distinct but related science investigations were awarded 
eight consecutive days (UT March 25 through April 1), totaling 
64 hr of clock time. However, problems associated with weather, 
observatory hardware, and high winds restricted our available 
time to the last two dates, resulting in a total of approximately 
10 hr of clock time on the comet. The lack of active guiding 
during daytime 7 also reduced the observing efficiency (e.g., 
see DiSanti et al. 2009 for details). Our measurement strategy 
balanced volatiles measured in all observing epochs (H 2 0, 
C 2 H 6 , and CH3OH) against those that were not measured with 
CRIRES or NIRSPEC (CO and CH 4 ), owing to inadequate 
Doppler shift on those dates. The resulting production rates 
and abundance ratios are listed in Table 2. 

NIRSPEC is a cross-dispersed echelle grating spectrometer 
that samples a wide range of spectral orders with relatively high 
spectral resolving power. On June 20, we used a 0?432 x 24" slit 
configuration (delivering X/ Ak ~ 25,000) for observations of the 


7 This resulted in an observing efficiency of approximately 37% for our 
CSHELL observations, based on 214 minutes of total on-source integration 
time over the two dates (UT March 31 and April 1; see Table 1). 
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Table 2 

Molecular Abundances for Volatiles in Comet C/2012 F 6 (Lemmon) a 


Species 

Date 

v b 

Lines c 

Trot 

Global <2 d 

Abundance 


(UT) 

(cm" 1 ) 


(K) 

( 10 27 s” 1 ) 

(Relative %) 

2013 Feb 2 e 
f « h = 1.20 AU; 4 

, = 0.99 AU; A-dot = 

: —2 km s _ 

1 ; P.A. = 262?4; a = 53? 4 




OH* 

Feb 2.31 

mult. 

24 

(70) 

187.88 ±27.56 

99 

CH 3 OH 

Feb 2.39 

2840.61 

10 

66 +1 9 lg 

3. 50 ± 0.47 

1.85 ±0.27 

c 2 H 6 

Feb 2.35 

2992.68 

6 

69+ 9 7 

0.49 ± 0.06 

0.26 ± 0.03 

c 2 h 2 

Feb 2.35 

3293.80 

9 

(70) 

< 0.10 

<0.05 

HCN 

Feb 2.31 

3299.52 

3 

(70) 

0.21 ± 0.08 

0.11 ±0.05 

h 2 co 

Feb 2.31 

2811.05 

7 

(70) 

< 0.22 

< 0.12 

nh 2 

Feb 2.31 

3310.10 

4 

(70) 

0.58 ±0.13 

0.31 ±0.07 

nh 3 

Feb 2.31 

3295.39 

1 

(70) 

< 1.20 

<0.63 

2013 Feb 4 
f R h = 1.18 AU; A 

= 0.99 AU; A-dot = 

—4 km s _1 

; P.A. = 287?6; a = 52?3 




h 2 o 

Feb 4.23 

3014.67 

8 

69+5 

189.05 ±11.04 

100 

OH* 

Feb 4.23 

3397.83 

3 

(70) g 

168.39 ± 14.98 

89 

HDO 

Feb 4.25 

2683.36 

11 

(70) 

<0.93 ±0.31 

<0.49 

2013 Mar 31 h 
f R h = 0.75 AU; A 

= 1.55 AU; A-dot = 

+16 km s _ 

1 ; P.A. = 199?3; a = 32?4 




h 2 o 

Mar 31.86 

mult. 

6/15 

80+\o 

459 ± 65 

100 

CO 

Mar 31.82 

2152.64 

2 

(80) 

19.5 ± 2.2 

4.25 ± 0.64 

c 2 H 6 

Mar 31.85 

2981.86 

2 

(80) 

1.80 ±0.24 

0.39 ± 0.07 

ch 3 oh 

Mar 31.85 

2984.34 

1 

(80) 

6.16 ± 0.68 

1.34 ±0.24 

CH 4 

Mar 31.93 

3040.64 

1 

(80) 

3.06 ± 0.57 

0.67 ±0.11 

h 2 co 

Mar 31.98 

2783.60 

3 

(80) 

2.46 ± 0.45 

0.54 ±0.12 

2013 Apr l h 
f R h = 0.75 AU; A 

= 1.56 AU; A-dot = 

+15 km s _ 

1 ; P.A. = 201?8; a = 31?8 




h 2 o 

Apr 1.90 

2984.34 

6/8 

100-i°o 

661 ± 92 

100 

CO 

Apr 2.03 

2152.64 

2 

(100) 

25.2 ± 2.3 

3.82 ± 0.63 

c 2 H 6 

Apr 1.81 

2985.00 

2 

(100) 

2.18 ± 0.22 

0.33 ±0.13 

ch 3 oh 

Apr 1.84 

mult. e 

3 

(100) 

6.76 ±0.51 

1.02 ±0.39 

2013 Jun 20 
{ R h = 1.74 AU; A 

= 1.79 AU; A-dot = 

+5 km s -1 

; P.A. = 256?3; a = 33? 4 




h 2 o 

Jun 20.50 

3441.46 

11 

U\ 

0 

1 + 
3V "- 1 

105.4 ± 6.6 

100 

c 2 h 6 

Jun 20.50 

2945.67 

13 

59 +14 

Jy -10 

0.31 ±0.03 

0.29 ± 0.03 

ch 3 oh 

Jun 20.50 

2908.89 

17 

50+5 

1.56 ±0.14 

1.48 ±0.16 

HCN 

Jun 20.55 

3299.54 

5 

(50) 

0.23 ± 0.02 

0.22 ± 0.03 

nh 3 

Jun 20.55 

3295.42 

1 

(50) 

0.64 ± 0.22 

0.61 ± 0.21 

nh 2 

Jun 20.55 

3301.84 

2 

(50) 

0.12 ±0.03 

0.12 ±0.03 

c 2 h 2 

Jun 20.55 

3295.94 

8 

(50) 

<0.06 

<0.06 

h 2 co 

Jun 20.55 

2851.41 

6 

(50) 

<0.25 

<0.24 


Notes. 

a Uncertainties represent lcr, and upper limits represent 3cr. The reported error in production rate includes the line-by-line scatter in 
measured column densities, along with photon noise, systematic uncertainty in the removal of the cometary continuum, and (minor) 
uncertainty in rotational temperature. 

b Mean wavenumber of all emission lines (used for this reduction) from a particular species. 

c “Lines” refers to the number of spectral intervals sampled. For the CSHELL water observations, the first entry indicates the number 
of intervals in the H20_3A setting that were used to determine T rot . The second entry indicates the total number of intervals used 
to calculate (9(H 2 0), including water lines in the H20_3A and CO_K settings, and OH* in the CH4_E and H2CO_B settings on 
March 31 and in the CH30H_A setting on April 1, as indicated in Table 1. 

d Global production rate, after applying a measured growth factor of (1.7 ± 0.2) and (1.7 ± 0.1) to the nucleus-centered production rates 
of volatiles measured on UT 2013 February 2 and February 4, respectively. For CSHELL (March 31 and April 1), growth factors were: 
H 2 0 (2.0 ± 0.2), CO (1.9 ± 0.1), C 2 H 6 (1.8 ± 0.05), CH 3 OH (2.3 ± 0.1), CH 4 (2.2 ± 0.1), andH 2 CO (1.8 ± 0.2). For NIRSPEC (June 
20), GFs were: 3.1 ± 0.2 for H 2 0 and 1.8 ±0.1 for minor species (based on measurements of C 2 H 6 , CH 3 OH, and HCN). 
e We adopt Q( H 2 0) from February 4, which is consistent with production rates from (mean) Q{ H 2 0) and (7(OH*) from the combined 
CH, HCN, and CH 3 OH settings on February 2. 

f R^-. heliocentric distance; A: geocentric distance; A-dot: geocentric velocity; P.A.: position angle of the extended Sun-comet vector; or. 
solar phase (Sun-comet-Earth) angle for the midpoint of data acquisition. 

g We adopted Trot (shown in parentheses) based on values measured for C 2 H 6 , CH 3 OH, and H 2 0 from CRIRES and NIRSPEC spectra, 
and for H 2 0 from CSHELL spectra when calculating the NC production rates for species whose T rot determination was not possible. 
h For trace species measured on both dates, we obtained the following weighted mean values, 10 27 s -1 (%): C 2 H 6 : 2.0 ± 0.2 (0.38 ± 
0.06), CH 3 OH: 6.0 ± 0.4 (1.25 ± 0.20), CO: 2.2 ± 0.2 (4.03 ± 0.45). 
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Table 3 

Summary of Molecular Abundances Ratios in Comet C/2012 F6 (Lemmon) 


Species 



Date (UT 2013) 


Weighted Mean 
% (H 2 0 = 100) 

Feb 2 

Feb 4 

Mar 31 

Apr 1 

Jun 20 

c 2 h 6 

0.26 ± 0.03 


0.39 ± 0.07 

0.33 ±0.13 

0.29 ± 0.03 

0.29 ± 0.02 

CH 3 OH 

1.85 ±0.27 


1.34 ±0.24 

1.02 ±0.39 

1.48 ±0.16 

1.48 ±0.11 

HCN 

0.11 ±0.05 




0.22 ± 0.03 

0.19 ±0.03 

CO 



4.25 ± 0.64 

3.82 ± 0.63 


4.03 ± 0.45 

CH4 



0.67 ±0.11 



0.67 ±0.11 

h 2 co 

<0.12 


0.54 ±0.12 


<0.24 

0.54 ±0.12 

nh 3 

<0.63 




0.61 ± 0.21 

0.61 ±0.21 

nh 2 

0.31 ± 0.07 




0.12 ±0.03 

0.15 ±0.03 

c 2 h 2 

<0.05 




<0.06 

<0.05 

HDO 


<0.49 




<0.49 


comet, with the slit positioned along the projected Sun-comet 
direction. Flux calibration was achieved through observations 
of the IR standard star HR 6707 with a (L720 x 24" slit 
configuration. 

Comet spectra were acquired in our standard four- step se- 
quence (ABBA), with accumulated times of 30, 60, or 120 s 
per stored A or B frame (depending on the instrument used and 
wavelengths targeted). In both the CRIRES and CSHELL cases, 
the telescope was nodded along the slit by 15" between the A 
and B positions, and with NIRSPEC the telescope was nodded 
by 12". 

Data reduction and analysis of the spectral frames included 
flat fielding, removal of pixels affected by high dark current 
and/or cosmic-ray hits, spatial and spectral rectification, and 
spatial registration of individual A and B beams. After combin- 
ing the A and B beams from the (processed) difference frames, 
which further assists in removing any residual background emis- 
sion, we extracted spectra by summing 15 spatial pixels (1"3) 
centered on the nucleus for the CRIRES observations, 10 pixels 
(2"0) for the CSHELL observations, and 9 pixels (1"3) for the 
NIRSPEC observations. We synthesized a transmittance func- 
tion for the terrestrial atmosphere by fitting the absorptions ob- 
served in the standard star spectra, and adjusting the estimated 
atmospheric contribution until an optimum fit was obtained. 
We used a multiple layer atmosphere and the LBLRTM model 
(Clough et al. 2005) that accessed the HITRAN 2008 molecular 
database modified with our custom updates (see, e.g., Villanueva 
et al. 2011b for details). We convolved our synthetic transmit- 
tance function to the resolving power of the comet observations, 
scaled it to the cometary continuum, and subtracted it from the 
measured spectrum. This isolated the net emission spectrum in 
excess of the cometary continuum. We then generated a fully 
resolved transmittance function for use in extracting the true 
transmittance experienced by individual spectral lines at their 
Doppler- shifted wavelengths. 

These observations allowed robust measures of pre-perihelion 
rotational temperatures and production rates for multiple species 
(Paganini et al. 2013b). The simultaneous measurement of H 2 0 
and its direct proxy (OH*) in each instrument setting permitted 
accurate characterization of the abundance ratio for each trace 
species. Spectra of H 2 0 and OH* are shown in Figure 1, 
and of other primary volatiles in Figure 2. Spectra resulting 
from the HDO setting are shown in Figure 3, along with a 
synthetic spectrum of expected HDO lines. With CRIRES, six 
volatiles (H 2 0, OH*, C 2 H 6 , HCN, CH 3 OH, and NH 2 ) were 
securely detected, and upper limits were retrieved for H 2 CO, 
NH 3 , C 2 H 2 , and HDO. Three of these species (H 2 0, C 2 H 6 , and 
CH 3 OH) were also measured with CSHELL and NIRSPEC. 


With CSHELL we also detected H 2 CO, CO, and CH 4 on March 
31/April 1 (Figure 4; DiSanti et al. 2013b). The latter two 
molecules require a geocentric velocity (A-dot) of ~10 km s -1 
or higher to displace comet lines from their (opaque) terrestrial 
counterparts (Figure 4). Comet C/2012 F6 presented small 
Doppler shifts during the CRIRES and NIRSPEC observations, 
preventing detection of CO and CH 4 (Table 1). 

3. RESULTS 

We determined the rotational temperature (T rot ) for an indi- 
vidual species, and we use this to quantify the production rate 
from the measured line fluxes. Results are given in Table 2, 
and the abundances obtained during our entire observing cam- 
paign are summarized in Table 3. The detailed procedures are 
discussed below. 

3.1. Rotational Temperatures 

The rotational temperature is obtained for an individual 
species by comparing modeled and observed fluxes for lines 
that span a sufficient range of rotational energies. The modeled 
spectral line intensities are derived from our custom quantum 
mechanical fluorescence models, which provide fluorescence 
efficiencies (g-factors) at the specified rotational temperature. 
The test temperature is modified iteratively until a satisfactory fit 
is obtained, and that optimum value for T mi (and its uncertainty) 
then permits a robust measure of production rate. The number 
of lines used for each species and the values retrieved for T mi 
and its uncertainty are as follows. With CRIRES, we measured 
r rot = 69+y K for ethane (6 emission lines) and 66 +1 9 ' K for 
methanol (10 lines) on February 2, and 69+ 8 5 K for H 2 0 (8 lines) 
on February 4. With CSHELL, we obtained T rot = 80 d= 10 K 
on March 31 and 100 zb 10 K on April 1, based on a setting 
containing six lines of H 2 0. Observations with NIRSPEC on 
June 20 resulted in T mt = 50 + J 6 K, 59+\ 4 0 K, and 50+ 7 K for 
water (11 lines), ethane (13 lines), and methanol (17 lines), 
respectively. See also Table 2. 

3.2. Production Rates and Molecular Abundances 

Using the optimum value for T rot , an apparent production 
rate is determined from the flux of each rovibrational tran- 
sition detected within our sampled aperture, Fy me (W m -2 ). 
Our methodology considers the molecular lifetime r(s), the 
molecular fluorescence efficiency or g-factor (g t ) at the ap- 
propriate Tj-ot and heliocentric velocity (u), terrestrial trans- 
mittance (7)), geocentric distance (A), and the fraction of the 
total molecular content in the coma sampled by each pixel f(x). 
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Comet C/2012 F6 (Lemmon) 

UT 2013 February 4 


H 2 0-ortho (t) 
M 2 U H 2 0-para (§) 



Wavenumber (cm 1 ) 

Figure 1 . Water and OH prompt emission (OH*) lines detected in comet C/2012 F6 (Lemmon) on UT 2013 February 4. H 2 O ortho and para emission lines are 
identified in the spectra. Here and in Figures 2 and 3, the gray region identifies the ±lcr uncertainty, and the underlying colored lines depict the modeled spectra. 

(A color version of this figure is available in the online journal.) 


The nucleus-centered production rate is calculated using 
„ _ 4jrA 2 fii ne 

— rji _C / \ * 

T 8iTj f(x) 

Our IR observations are impacted by slit losses (i.e., loss of 
flux) that result from atmospheric “seeing” and slight aperture 
effects (e.g., imperfect centering of the cometary photocenter in 
the slit). To correct for these effects, we obtain a growth factor 
(GF) from the observed spatial profiles after averaging Q at 
diametrically opposite positions along the slit, as this averages 
asymmetries in outflow (Xie & Mumma 1996). For our infrared 
observations, multiplying the Q nc by GF results in total (or 
“global”) production rates, Q tot . 

Details of fluorescence models for each molecule can be 
found as follows: H2O, HDO (Villanueva et al. 2012b), OH* 
(Bonev et al. 2006), C 2 H 6 (Villanueva et al. 2011a), CO, C 2 H 2 , 
and CH 4 (Paganini et al. 2013a; Villanueva et al. 2011b; Gibb 


et al. 2003), NH 3 , HCN (Villanueva et al. 2013; Lippi et al. 
2013), NH 2 (Kawakita & Mumma 2011), H 2 CO (DiSanti et al. 
2006), and CH3OH (Villanueva et al. 2012a; DiSanti et al. 
2013a). 

Obtaining reliable spatial profiles for observed volatiles re- 
quires adequate signal-to-noise ratios. CRIRES achieved such 
profiles for H 2 0 and C 2 H 6 , and for the continuum (see 
Figure 5(A)). These spatial profiles (pre-perihelion) reveal a 
fairly symmetric distribution of outgassing around the nucleus 
out to projected distances near 2000 km (as projected onto the 
sky plane), with some excess in the range ~ 1000-1500 km, per- 
haps indicating the presence of icy grains toward the northwest 
(i.e., relative anti-sunward direction). 

With CSHELL (Figure 5(B)), we obtained spatial profiles 
for H 2 0, C 2 H 6 , CO, and H 2 CO. Compared to profiles obtained 
with CRIRES, the spatial distribution of H 2 0, C 2 H 6 , CO, and 
the continuum extended to significantly larger nucleocentric 
distances (to ^4000 km), owing to the larger production rates at 
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Comet C/2012 F6 (Lemmon) 

UT 201 3 February 2 


1000 



Wavenumber (cm 1 ) 

Figure 2. CRIRES observations of minor volatiles on UT 2013 February 2. Panel (A) displays the methanol V 3 band and several OH* emission lines. Panels (B) and 
(C) show the 2982-3012 cm -1 range covering the ethane V 7 band (P, R, and Q lines) along with a line of OH* and several emission lines of methanol. Panel (D) 
shows HCN emission lines with radical NH 2 . 

(A color version of this figure is available in the online journal.) 


0.75 AU versus 1.2 AU (with CRIRES). Formaldehyde extended 
to -3500 km and displayed stronger intensities toward the anti- 
sunward direction. The continuum showed a rather symmetric 
distribution about the nucleus, while CO was enhanced in the 
sunward direction and C 2 H 6 showed some excess flux in the 
anti-sunward direction. In agreement with CRIRES, the H 2 0 
profile is rather symmetric with some enhancement in the anti- 
sunward direction at about 2000 km, similar to C 2 H 6 . The 
different configurations of CO, H 2 0, and C 2 Hg could be related 
to separate polar and apolar release of ices in the nucleus (see 
Section 4.3). NIRSPEC profiles also revealed enhancement of 
C 2 H 6 in the anti-sunward direction (similar to HCN and H 2 0; 
Figure 5(C)). On the other hand, methanol shows a rather 


symmetric distribution about the nucleus, and the continuum 
was stronger in the sunward direction. The observed (flux) 
enhancement of H 2 0 relative to other minor volatiles could 
be related to the release of icy grains at = 1.78 AU. 

3.2.1. CRIRES Observations 

The production rates and mixing ratios for water and seven 
trace species are given in Table 2. Our value for water production 
on February 4 (g tot = 1.9 ± 0.1 xlO 29 molecules s -1 ) agrees 
with estimates of H 2 0 production using OH prompt emission 
lines (OH*, a direct proxy for water) detected in the CH, 
CH 3 OH, and HCN settings on February 2, and in the HDO 
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— HDO(xio) 

Comet C/201 2 F6 (Lemmon) Model for HDO at 70 K 



I I I I I I I I I I I I I I I I I I I I I 

2666 2664 2662 2660 2658 2656 




Figure 3. Spectra resulting from the CRIRES HDO setting. Observations of HDO were strongly affected by terrestrial extinction due to the low geocentric velocity 
of C/2012 F6 on February 4 (—2 km s -1 ), and this limited our sensitivity. 

(A color version of this figure is available in the online journal.) 


and H2O settings on February 4 (Table 2). Our results show 
no evidence for significant variation in water production on 
February 2 and 4; we therefore conclude that adopting this 
value for both nights provides an accurate measure of molecular 
abundances. 

3.2.2. C SHELL Observations 

CSHELL observations occurred shortly after perihelion, 
when the gas production was larger than when measured pre- 


perihelion with CRIRES. The water production rate on April 1 
(6.6 ±0.9 x 10 29 s- 1 ) was marginally higher than on March 3 1 
(4.6 zb 0.6 x 10 29 s _1 ), an increase by a factor of 1.5 ± 0.4, 
or ~ 1.5 o based on their combined uncertainties. The produc- 
tion rates for C2H6, CH3OH, and CO were also somewhat 
higher on April 1, but the abundance ratios for trace gases 
were systematically lower on April 1 than on March 31 (see 
Table 2). However, both quantities were unchanged within their 
lcr uncertainties, even though the slit was oriented north-south 
on March 31 and east-west on April 1 (see Table 2). This 
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C/2012 F6 (Lemmon) 

IRTF/CSHELL, UT 2013 March 31/April 1 2 3.3° Solar elongation 



2845 2844 2843 2842 2841 2840 2839 

Wavenumber (cm 1 ) 


Figure 4. Example of volatiles detected with CSHELL along with model fits of the cometary continuum (dashed lines) to extract emission spectra. Our fluorescence 
models of these spectra allow us to obtain production rates for these species. Panel (A) shows the 2150-2156 cm -1 range, in which the continuum- subtracted residuals 
contain the emission lines of H 2 O and CO. Similarly, panels (B) and (C) show detections of H 2 CO, CH 3 OH ( Q branch), and OH*. 

(A color version of this figure is available in the online journal.) 


suggests a non- varying composition for the responsible source 
region(s). 

3.2.3. NIRSPEC Observations 

Our observing campaign concluded with a run on June 20, 
when comet C/2012 F6 was at 1.74 AU from the Sun and 
1 .79 AU from Earth. The water production rate was 1.1 zb 0. 1 x 
10 29 s-\ and we securely detected C 2 H 6 , CH 3 OH, HCN, NH 3 , 
and NH 2 . The results are listed in Table 2. 


4. DISCUSSION 
4.1. Water Production 

The H 2 0 production rates we measured with CRIRES and 
CSHELL follow a heliocentric power law of R h “ 2 2±0 8 , which 
is consistent with insolation-limited outgassing, while a com- 
parison of CSHELL and NIRSPEC observations resulted in a 
power law of Rj/ 1 9±0 4 . A possible explanation for the relative 
increase in water production between March 3 1 and April 1 is 
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VLT/CRIRES 

A R h = -1.20 All 

UT 2-4 February 201 3 


IRTF/CSHELL 

B R h = 0.75 AU 

UT 31 March - 1 April 2013 


Keck-2/NIRSPEC 

C R h = 1 .74 AU 

UT 20 June 2013 



2000 1000 0 -1000 -2000 
Nucleocentric distance (km) 


2000 0 -2000 
Nucleocentric distance (km) 



Nucleocentric distance (km) 


Figure 5. One-dimensional spatial profiles of primary volatiles and the continuum during the inbound (A: CRIRES) and outbound (B: CSHELL, C: NIRPSEC) 
passage of comet C/2012 F6 (Lemmon). The cardinal directions (i.e., the Sun-comet vector and solar phase angle) are indicated relative to the slit orientation, which 
is horizontal with respect to these plots. For further details, see Section 3.2. 

(A color version of this figure is available in the online journal.) 


the presence of strong jets near perihelion; however, in this case 
we would expect pronounced asymmetries that are not seen in 
our H 2 0 profiles (Figure 5). The most straightforward expla- 
nation favors an increase in overall gas production on April 1, 
perhaps due to preferential exposure of a larger or more active 
vent as a result of nucleus rotation, although the limited data 
sample restricts any conclusive interpretation. 

In Figure 6, we compare our measured H 2 0 production rates 
with those inferred from hya using the SWAN all-sky camera 
on the Solar and Heliospheric Observatory (SOHO; Combi 
et al. 2013), and preliminary results inferred from the optical 
OH using the TRAPPIST telescope (revised estimates; Opitom 
et al. 2013a, 2013b; C. Opitom 2013, private communication). 
We do observe a systematic difference between the TRAPPIST 
results (from measurements of OH using 2(H 2 0) = 1.1 g(OH)) 
and those obtained with other techniques. 

We also compare with g(H 2 0) estimated from visual mag- 
nitudes from Seiichi Yoshida’s webpage 8 using the empirical 
relation log 2(H 2 0) = 30.675 — 0.2453(m v — 5 log A) given 
by Jorda et al. (2008). In early February, all three techniques 
(CRIRES, SWAN, and empirical estimates) provided similar re- 
sults. However, the Jorda relation should be used with caution, 
since it does not account for the large differences seen among 
comets in C 2 emission intensities, in grain size distributions, 
and also in dust/gas ratios, all of which can influence the visible 
magnitude and thus introduce significant systematic uncertain- 
ties. Bearing in mind these limitations, these optical magnitudes 
provide estimates of g(H 2 0) in the range of 1.2-2. 2 x 10 29 s -1 , 
consistent with SOHO / SWAN values (1. 6-2.1 x 10 29 s -1 ) and 
also with our CRIRES value (1.9 zb 0.1 x 10 29 s -1 ). 

Our CSHELL production rate for H 2 0 on March 31 (4.6 zb 
0.6 x 10 29 s- 1 ) is lower than the prediction based on the 
Jorda estimate (~(6.2-11.2)x 10 29 s -1 ); however, our value on 
April 1 (6.6 zb 0.9 x 10 29 s -1 ) is consistent with the predicted 
value at the 1 o level. The SWAN values show better agreement 
with the predicted estimates. However, although each SWAN 
measurement represents the mean water production averaged 
over ~ 10 days, the day-to-day variations near perihelion exhibit 
significant excursions — relatively larger than the (marginal) 
increase we measured between March 3 1 and April 1 . Similar 

8 http://www.aerith.net/comet/catalog/2012F6/2012F6.html 
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Figure 6. Comparison of measured and estimated H 2 0 production rates: IR 
facilities VLT/CRIRES, IRTF/CSHELL, and Keck2/NIRSPEC; the SWAN 
all-sky Ly a camera on the Solar and Heliospheric Observatory (Combi et al. 
2013); the TRAPPIST telescope (revised estimates; Opitom et al. 2013a, 2013b; 
C. Opitom 2013, private communication); and empirical estimates. The latter 
estimates represent only first-order approximations of <2(H 2 0) derived from 
visual magnitudes (m v ) at optical wavelengths and the empirical relation of 
Jorda et al. (2008), log <2(H 2 0) = 30.675 — 0.2453(m v — 5 log A). For m v we 
used Yoshida’s estimates based on optical observations(m v = 5.2 + 10 log Rh); 
see http://www.aerith.net/comet/catalog/2012F6/2012F6.html. We arbitrarily 
include an uncertainty crm w ± 0.5, displayed as the shaded gray region. These 
estimates are subject to certain uncertainties discussed in Section 4.1; however, 
the shape of the curve (variation with heliocentric distance) is less subject to 
such variations. 

(A color version of this figure is available in the online journal.) 
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variations are not seen in the SWAN results in early February, 
around the time of our CRIRES observations. This could point 
to substantially more variable gas production around perihelion. 

The water production rate we obtained with NIRSPEC (1 . 1 zb 
0.1 x 10 29 s -1 ) on June 20 (Rh = 1-74 AU) allows us to quantify 
a rather comprehensive evolution pre- and post-perihelion of 
comet C/2012 F6. At the 2 a level, our measurements display 
a behavior similar to estimates derived from the SOHO / SWAN 


9 


The Astronomical Journal, 147:15 (llpp), 2014 January 


Paganini et al. 


results and the Jorda approximation. Even though our IR 
observations provide a general overview of the overall evolution 
of water production and abundance of minor volatiles (relative 
to water), upcoming measurements from other facilities will add 
key information to the evolution of water production rates, and 
thus the real behavior of comet C/2012 F6. 

4.2. Abundances of Trace Volatiles — Placing 
C/ 2012 F6 ( Lemmon ) in Context 

Our water production rates on March 31 and April 1 were 
a factor of ~2-3 higher than on February 4, which is reason- 
able considering the difference in R h (Table 1). The observed 
decrease between March 31 /April 1 and June 20 is likewise 
consistent with increasing heliocentric distance. For a homo- 
geneous volatile composition, the abundance ratios of minor 
volatiles relative to H2O are expected to be independent of he- 
liocentric distance — at least within the regime of H 2 0-driven 
activity encompassing our observations. 

We observed only two minor volatiles (C2H6 and CH3OH) 
in all three observing runs (Tables 1-3). The abundance ratios 
for C 2 H 6 and CH 3 OH relative to H 2 0 on February 2 agree 
(within 1.5a) with their corresponding values on March 31 and 
April 1, and also (within la) with values on June 20, suggesting 
consistent relative production of these (key) species pre- and 
post-perihelion. 

Based on the weighted mean abundances (Table 3), our 
observations show that C2H 6 in C/2012 F6 is depleted by 
approximately a factor of two relative to that measured in the 
majority of OC comets (i.e., relative to the current “normal” 
value; see Mumma & Charnley 2011), and similarly for CH3OH 
(typical “normal” values for C2H6 and CH3OH are ~0.5% and 
~2%-3%, respectively). Our measured abundance of CO is 
consistent with its median value found among comets (~4%), 
HCN is consistent with its “normal” abundance (~0.2%), CH4 
is slightly below its median value (0.9%-l%), and H2CO is 
above the values found for (native) formaldehyde in a number 
of comets measured to date (typically ~0.1%-0.2%) and is 
closer to the higher abundances (>0.5%) seen in C/2002 
T7 (FINEAR), C/2002 Cl (Ikeya-Zhang), and 9P/Tempel 1 
(DiSanti et al. 2002, 2006; Mumma et al. 2005). 

4.3. Possible Interpretation of Chemical 
Abundances in C/ 2012 F6 

The observed depletion of some volatiles in this comet could 
indicate formation in, or a possible (later) exposure to, warm 
environments, whose temperatures were high enough to produce 
subsequent redistribution of these volatiles within (or escape 
from) the nucleus. However, the normal abundance of HCN, 
CH 4 , and CO, plus the enhancement of formaldehyde (believed 
to have formed beyond the CO snow line via H-atom addition to 
CO-rich ices; Tielens & Hagen 1982), argues against formation 
under such conditions and suggests possible heterogeneity in 
this comet. The latter suggests that either (1) there was radial 
mixing during the formation of this icy body, (2) some volatiles 
were selectively shielded in water ice mantles during accretion, 
or (3) that the nucleus is composed of different structural phases 
of ice (polar and apolar), as, for instance, observed recently in 
comets C/2007 W1 (Boattini), 103P/Hartley 2, and C/2009 PI 
(Garradd; Villanueva et al. 2011a; Mumma et al. 2011; Paganini 
et al. 2012). Alternatively, post-processing in the protoplanetary 
disk (cf. the Grand Tack model) or in the cometary reservoirs 
(OC/KB) should not be disregarded. 


4.4. D/Hin Water 

Observations of HDO were adversely affected by terrestrial 
extinction due to the low geocentric velocity of C/2012 F6 
on February 4 (—2 km s -1 ), and this limited our sensitivity 
(Figure 3). Regardless of these limitations, we obtained an upper 
limit (3a) for HDO of <0.9 x 10 27 s -1 , corresponding to a 
D/H ratio in water of <2.45 x 10 -3 (< 16 VSMOW). Although 
not restrictive, this upper limit is consistent with values found in 
other comets (Mumma & Charnley 2011 and references therein; 
Hartogh et al. 2011; Bockelee-Morvan et al. 2012; Fis et al. 
2013). 

4.5. Ortho-Para Ratio in Water 

If we assume that the ratio of nuclear spin species remains 
unaltered after incorporation in the cometary nucleus, the 
ortho-para ratio (OPR) of water should serve as an indicator 
of thermal conditions during its origins. A measurement of the 
OPR was obtained from water emission lines detected using 
CRIRES observations on February 2 and with NIRSPEC on 
June 20. Figure 1 shows the CRIRES detection of ortho- and 
para- water and identifies these lines along with OH* . To estimate 
the OPR in comet C/2012 F6, we used eight ortho- water and 
seven para- water emission lines on February 4 (with CRIRES) 
and six ortho- water and four para- water emission lines on June 
20 (with NIRSPEC). The resulting OPRs (3.78 ± 0.77 with 
CRIRES and 3.20 ± 0.37 with NIRSPEC) agree within la. 
Their weighted mean OPR is 3.31 ±0.33, which constrains the 
spin temperature to be >37 K at the 95% confidence limit but 
is consistent with statistical equilibrium. For further details on 
the significance of OPR measurements, we refer the reader to 
Mumma & Charnley (201 1), Bonev et al. (2013), and references 
therein. 

5. SUMMARY 

We investigated the production of primary volatiles in comet 
C/2012 F6 (Femmon) using VFT/CRIRES on 2013 February 2 
and 4 (pre-perihelion; R h = 1.2 AU), IRTF/CSHEFF on March 
31 and April 1 (post-perihelion; R h = 0.75 AU), and Keck- 
2/NIRSPEC on June 20 (post-perihelion; R h = 1.74 AU). We 
detected 10 volatile species (H 2 0, OH* prompt emission, C 2 H 6 , 
CH3OH, H 2 CO, HCN, CO, CH 4 , NH 3 , and NH 2 ), and obtained 
upper limits for two others (C 2 H 2 and HDO). We quantified 
the OPR of H 2 0 using ortho and para emission lines (OPR = 
3.31±0.33), which constrains the spin temperature to be >37 K 
at the 95% confidence limit but is consistent with statistical 
equilibrium. Our measured production rate for H 2 0 and upper 
limit for HDO provided an upper limit for D/H in water 
(D/H < 2.45 x 10- 3 , i.e., <16 VSMOW). 

With CSHEFF we obtained a marginally higher water produc- 
tion rate on April 1 compared with March 3 1 , but the values on 
both dates showed a reasonable excess compared with CRIRES 
results in February (and with NIRSPEC results in June), in ac- 
cord with the increased (decreased) insolation between these 
epochs. 

The obtained volatile abundance ratios are summarized in 
Table 3. Comet C/2012 F6 is rather depleted in C 2 H 6 and 
CH3OH; HCN, CH 4 , and CO displayed normal abundances, and 
H 2 CO was the only volatile showing an enhancement relative 
to water. These volatile abundances suggest a certain degree 
of chemical post-processing and/or formation within warm 
regions in the nebula. However, the normal abundance of HCN 
and of the hypervolatiles CH 4 and CO, plus the enhancement 
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of H 2 CO, argues against formation under such conditions, and 
suggests possible heterogeneity in this comet. Further evidence 
on this matter will be revealed as additional results from radio, 
infrared, and optical wavelengths become available. 
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